A semiquantitative approach to the impurity-band-related transport properties of 
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We investigate the spin-polarized transport of GaMnAs nanolayers in which a ferromagnetic order 
exists below a certain transition temperature. Our calculation for the self-averaged resistivity takes 
into account the existence of an impurity band determining the extended ("metallic" transport) or 
localized (hopping by thermal excitation) nature of the states at and near the Fermi level. Magnetic 
| order and resistivity are inter-related due to the influence of the spin polarization of the impurity 

band and the effect of the Zeeman splitting on the mobility edge. We obtain, for a given range 
of Mn concentration and carrier density, a "metallic" behavior in which the transport by extended 
Q ' carriers dominates at low temperature, and is dominated by the thermally excited localized carriers 

near and above the transition temperature. This gives rise to a conspicuous hump of the resistivity 
which has been experimentally observed and brings light onto the relationship between transport 
and magnetic properties of this material. 



I. INTRODUCTION 



Diluted magnetic semiconductors (DMS) are important materials in Spintronics since, in their ferromagnetic phase, 
'"O ' they can be used as a source or a filter of spin-polarized carriers to be injected into semiconductor heterostructures. 
It is interesting to have the transition temperature as high as possible, aiming operation at room temperature. In 
the study of the possible origins for the magnetic order, the existence of an impurity band has been questioned. An 
impurity band occur when impurities, donors or acceptors, are close enough to have a finite width of the density 
of states (DOS) for the bound electron or hole. It is well known that at high concentration the maximum of DOS 
approaches the conduction band (donor impurities) or the valence band (acceptor impurities), and at the same time 
a tail appears in the extended state band (conductance or valence). At impurity concentrations still higher the two 
DOS coalesce, and the energy gap between then disappears. The Mott- Anderson transition occurs when the states at 
the Fermi level change their character from localized to extended. The Fermi level is said, then, to have crossed an 
edge of mobility, separating localized and extended states. Although, at least up to present, this transition cannot be 
classified as a phase transition with known critical exponents etc, the concept of the mobility edge is well accepted. 
If it plays an important role in the magnetic properties of DMS, it is a question of identifying somehow the existence 
of an impurity band itself. 

Gai-^Mn^As is a semiconductor alloy belonging to the DMS family. Example of DMS materials are among oxides 
(like ZnO doped with Eu, Co etc), II- VI compounds (like CdTe doped with Mn), III-V compounds (like GaMnAs), 
and IV- VI compounds (like PbTe doped with Eu) . The possibility of controlling the energy gap of the non-magnetic 
alloys GaAlAs and InGaAs, the material's n or p character by doping, the electron-hole recombination time, the 
effective mass, among other properties, together with the capability of growing these materials in heterostructures 
of high quality, makes GaMnAs special. After the work of van Esch et air- an d Oiwa et al^- where they produced 
good homogeneous samples by low-temperature molecular beam epitaxy and obtained ferromagnetic phases with 
high transition temperatures (Tc), many improvements brought this temperature to 180K in epilayers£ Lower Tc 
were obtained in multilayer system, which are important for devices.— Recently, much higher Tc near 250 K were 
obtained in highly doped digital layers^ Although higher transition temperatures have been obtained in other DMS, 
the characteristics given above for GaAs focus the attention on GaMnAs. Besides of its importance for Spintronics 
applications, it shows challenging questions concerning the understanding of its magnetic and transport properties. 

Taking as a basis for discussion those samples in Refs. [l] and d, a non-metal-to-metal transition is observed as 
the Mn concentration increases and goes above x — 0.03. This is followed by a metal-to-non-metal transition near 
x = 0.065. On the other hand, the ferromagnetic phase is observed even before the first transition, but in this 
case its ferromagnetic-to-paramagnetic transition temperature is low. At first this temperature increases with the 
Mn concentration, reaching a maximum inside the concentration range corresponding to the metallic phase. For 
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concentrations above the one giving rise to the metal-to-non-metal transition, ferromagnetism persists but with 
decreasing ferromagnetic-to-paramagnetic transition temperature. Recently, the quality of the samples increased still 
further, leading to the possibility of doping GaAs with higher concentrations of Mn and obtaining higher Curie 
temperatures. The conductive behavior, concerning the non-metal-to-metal transition and vice versa, together with 
its influence over the transition temperature are still observed in these samples. 3 *^ In Mn doped InAs quantum wells 
resistivity measurements have also shown localization effects in the hole system^ 

During the last ten years a lot of attention was given to understand the mechanism leading to the ferromagnetic order 
in GaMnAs. Several possibilities were explored, such as RKKY interaction^ mean field approximation, 10 magneto- 
polaronii etc. We have performed^ 2 " self-consistent Monte Carlo simulations using an RKKY-like model to obtain the 
magnetization, susceptibility and Curie temperature of several GaMnAs/GaAs heterostructures. We concluded that 
the possibility of creating an inhomogeneous spin-polarization density across the sample, more specifically, having 
a highly concentrated polarization density inside the magnetic layer, increases remarkably the Curie temperature. 
More recently, Green's function formalisms in different approaches were used to treat this problem. One of them 
uses a Matsubara Green's function and modeled the spin correlation by a Brillouin function to obtain the conductive 
behavior on ferromagnetic quantum dots.— Also, the Kondo lattice model was used to treat the diluted regime in 
the ferromagnetic phase in the DMS^ as well as a local random phase approximation^ All these methods consists 
of treating disorder in the magnetic interaction, but neglect Coulomb scattering by impurities. As a consequence, 
localization is not considered. On the other hand, a survey on the optical and transport properties of GaMnAs samples 
over a wide range of concentration s 16 ' 17 identified the existence of an impurity band and explored the character of 
the states near the Fermi levels This latter, based on photoluminescence measurements, found that even in samples 
of Gai-^Mn^As doped with very low Mn concentration, i.e., in non-metallic samples, there are evidences of the 
formation of an impurity band. For metallic samples with a doping > 2%, no effect of an isolated impurity band 
was perceived in terms of dc transport, suggesting a fusion of the impurity band with the valence band. Impurity 
bands in GaMnAs have been also suggested after weak localization observation^ Calculations of impurity band have 
also been performed i 20 ' 21 It is known that the character of holes determines the indirect interactio n 14 ' 22 among the 
magnetic ions, influencing the Curie temperature of the magnetically ordered samples. Therefore, our motivation 
in this paper is to explore two aspects that we consider as key issues to understand the inter-related magnetic and 
transport properties in GaMnAs: spin polarization and localization. 

For decades the nature of the states of the non-interacting and the interacting quasi-two-dimensional electron 
systems involving, for instance, the metal-to-non-metal transition, the possibility of having extended states in absence 
of an external magnetic field, and the existence of a minimum metallic conductivity crossed periods where things 
seemed to be settled, alternating with periods where new questions and observations appeared. For the sake of having 
a historical view of those questions, one can follow the discussions appearing in Refs. [23| and [24]. In the present case, 
carriers are holes, their concentration is much higher than in typical Mott- Anderson transition in semiconductors (in 
the order of 10 18_19 cm~ 3 , while in Gao.95Mno.05As it is around 10 20 cm~ 3 ). Besides, the spin-polarization in GaMnAs 
is due to the local potential interaction, generally accepted to be modeled by an sp-d Kondo-like potential, instead of 
being due to an external magnetic field. Disorder enters in the random Coulomb scattering by impurities, and also 
on the sp-d term, in the latter due to the stochastic interaction with the fluctuating local dipole moments together 
with the random location of these dipoles. 

Recently the effects of the disorder in a quasi-two-dimensional GaMnAs layer was treated using a self-consistent 
multiple-scattering approximation to calculate its density of states and spectral density function. Using parameters 
corresponding to GaMnAs thin layers, a wide range of Mn concentrations and hole densities were explored to unders- 
tand the nature, localized or extended, of the spin-polarized holes at the Fermi level for several values of the average 
magnetization of the Mn system. For a certain interval of Mn and hole densities, an increase on the magnetic order 
of the Mn ions comes together with a change of the nature of the states at the Fermi level, showing a derealization 
of spin-polarized extended states anti-aligned to the average Mn magnetization, and a higher spin-polarization of the 
hole gas. This was associated to the appearance of the non-metal-to-metal transition caused by the increase of Mn 
concentration, once the ferromagnetic phase is reached, and was observed experimentally since the first samples were 
produced. The metal-to- non-metal transition observed in samples of Refs. [l] and [U appears in that calculation as 
Anderson transitions occurring in the spin anti-aligned impurity band of samples with higher Mn concentration. 

In this paper we link the magnetic and the transport properties coming out of a picture in which the existence 
of an impurity band of spin-polarized holes determines the localized or extended nature of states at and near the 
Fermi level. At a given average magnetization, for a specific sample, the Zeeman-like splitting, as described below, 
determines a mobility edge. Then, two transport mechanism take place. One, a thermally activated hopping of the 
localized states, those occupied states lying below the mobility edge in the impurity band. Another, the "metallic" 
transport, corresponding to the contribution of the occupied extended states, whose energies lie between the mobility 
edge and the Fermi level. We have, therefore, two channels like parallel circuits, and if we describe the transport by 
the resistance, it will be given by law of the inverses. The calculation herein proceeds as follows. Section II shows the 
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method used to obtain the density of states, how we proceed to the classification of the states as extended or localized, 
and in consequence, how we obtain a mobility edge. Section III describes a method allowing directly the calculation of 
the self-averaged "metallic" resistance. Then, Section IV applies these formalisms to a double layer sample providing 
a highly inhomogeneous spin polarization density, seeking for an ideal combination of Mn concentration and density 
of holes. Finally Section V compares our results with experimental observations and makes general comments based 
on these results. 



II. QUASI-TWO-DIMENSIONAL IMPURITY BAND AND THE CHARACTER OF THE HOLE STATES 

Our system is described by the following Hamiltonian: 

H = Y^ e P a PCT a pCT + ^imp(q) a i+w a P° + *W(q, 0", a') ^ a P+q(T a p(T , . (1) 

pa qcr p qerer' p 

The second term on the right represents the Coulomb scattering by impurities, and the third term is the interaction 
of a hole with the localized magnetic moments at the Mn positions. It is modeled by a Kondo-like potential: 

V sp - d (f) =-lJ2Si- s(r)«5(r - % (2) 

i 

where the localized spin of the Mn ion Si at position Ri is treated as a classical variable, and s(r) is the spin operator 
of the carrier at position r; / is the sp~ d interaction. The third term on the right represents the Coulomb scattering 
by impurities, eventually the same ions giving rise to the local dipole moments. 

The impurity band is obtained by averaging the Green's function (GF) on the dipole moment orientation and on 
the impurity position, which, for a given scattcrcr configuration is defined by: 

G R (p, t; p', t') = -i (0|7> p (<)^, (t')]|0) . (3) 

Taking a diagrammatic expansion we reach a Dyson equation for the averaged GF by including all irreducible insertions 
to form the self-energy S i? (p,E). The configurational averaged GF (which is diagonal in p, as a consequence of the 
space becoming homogeneous after the averaging process) can be written in terms of the Dyson equation: 

< G*(p; E) >= G£(p; E) + G«(p; E)E*(p, E) < G*(p; E) > . (4) 

Obtaining this average is an extremely complicate task, which has not been achieved yet. Our choice has been to treat 
separately the magnetic and Coulomb scattering. The first one is taken in zero-th order, neglecting fluctuations of the 
magnetic moments on the impurities, which otherwise produce spin-flip scattering. Besides, we neglect fluctuation on 
the magnetization and assume homogeneous distribution of the localized magnetic moments with a density given by the 
Mn concentration with an average magnetization < S >— 5H/2 < M >, where (M) is the normalized magnetization 
(0 < (M) < 1) . This leads to a real contribution to the self-energy which is p- independent, and can be incorporated 
to the energy particle 

Up) = co(p) " ^N /3(M)a, (5) 

This approximation corresponds to the zero-th order term in Ref. fl3l No fix = Ic, c stands for the Mn density, 
a = ±1, x is the Mn doping factor, and Nq(3 is the exchange potential for holes, Nq(3 — —1.2eV, according to Ref. 

Next we proceed to the calculation of the self-energy due to impurity scattering. In order to simplify the notation 
we will omit the subscript a, since the derivation is performed for each spin orientation individually. Due to the high 
impurity density, the calculation of the average Green's function must consider terms beyond the self-consistent Born 
approximation. This problem has been handled successfully in the past with the use of the multiple-scattering appro- 
ximation which consists in considering self-energy insertions due to repeated scattering by the same impurity. Here 
we use the method developed by Klauder— with the improvements performed by Ghazali and Serre.— The multiple 
scattering approximation consists in selecting from the self-energies only those terms described by the diagrammatic 
expansion in Fig. [TJ expressed by the expansion: 

S el (p; E) = -^w / Ai<qi)G°(p + qi)u(-qi) 

+ (2^p / ddqi J ddc l2«(qi)G (p + qi Mq 2 - qi)G°(p + q 2 )w(-q 2 ) + ... (6) 
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Next, we define the vertex function K(k, qi; E), which obeys the following Dyson equation: 

K(k, qi ;E) = J (f'qi^cu -qi)G^)[2V«(k-qi) + A-(k,qi;£)]. (7) 

It can be shown, as derived in Refs i 27 i 28 , that the self-energy can be obtained from the vertex function by: 

E ei (k;£) = K(k,k;E). (8) 

So, S e i(k; E) is the diagonal term corresponding to line k, row k of matrix K, calculated at energy E. For the sake 
of making the formalism operational, a function U is defined such that: 

U(k,q;E) = K(k,q;E)+Nv(k-q). (9) 

In that case, Eq.flT]) becomes: 

U(k, q; E) = Nv(k - q) + J d d qWq' ~ q)G(^)U(k, q'; E). (10) 

With this result we have a linearized matrix equation: 

[I - vG]U = Nv, (11) 

where the sign tilde is used to identify a matrix. Therefore, the problem is reduced to a simple linearized matrix 
equation of the kind A ■ X — B, where B — Nv is a fixed matrix, while A — [I — Gv), and X = U changes at each 
iteration in the self-consistent process. 

If x = 0.05, Vmag introduces a kind of Zeeman splitting of 150 meV for fully magnetic ordered samples. It is 
important to emphasize the fact that the splitting in two bands, one with spins aligned with the average magnetization 
(which goes up in energy) and the other with spins anti- aligned (which goes down in energy), is not really a Zeeman 
splitting since it does not result of the interaction of the spin with an external magnetic field, but it comes out of a 
configurational average performed on the sp-d interaction of the Kondo-like local potential. 

The self-consistent calculation is performed for a two- dimensional hole gas of areal density n s submitted to Coulomb 
scattering by a negative ionized impurity system of concentration ATj. These two data, n s and iVj, are considered 
as independent parameters. This is important in the present context, since it is known that the density of free 
carriers in the ferromagnetic GaMnAs samples is just a fraction of the concentration of Mn, the reason for this being, 
still, a controversial question in the literature. Besides, we leave the method open to the possibility of doping with 
non-magnetic dopants, or changing the carrier concentration by applying a gate voltage. We start the self-consistent 
calculation with a free-particle spin-polarized Green's function G° (k, E) : 

* e ^-,i) ti .> ' (12) 

Next, we calculate the density of states as a function of the impurity and carrier concentrations. A deeper unders- 
tanding on the localized/extended character of the states can be attained by analyzing the spectral density function 
(SDF), 

A a (k, E) = --Im[< G*(k; E) >] (13) 

at a given energy as a function of the wave vector. Then, we observe the typical shapes corresponding to localized 
and extended states, as described in the following discussion. 

An electron bound to an impurity in the very diluted regime has a bell-shaped SDF with the maximum at k = 0. 
A free particle, having its energy perfectly determined by the wave vector, E(k) — ti 2 k 2 /2m has its SDF at a given E 
represented by delta function in variable k centered at k = y2m*E/H. When random scatterers are included, a width 
is introduced in the density of states of the single impurity bound state, together with a tail in the bottom (top) of the 
conduction (valence) band. As the impurity concentration increases, the isolated impurity band merge into the tail, 
and the gap between impurity states and free carriers disappear. The shape of the SDF for a given energy is used as a 
semi-quantitative criteria for the identification of a state at energy E as being extended or localized. For instance, if 
the energy is deep inside the conduction/ valence band the SDF plotted as a function of the wave vector k shows a very 
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sharp distribution at a given k, while a state deeply inside an impurity band appears as a broad Lorentzian centered 
at k = 0. This is schematically shown in Fig. [5J At left, (a) shows the DOS of a characteristic two-dimensional 
n-doped semiconductor system with n s = 2 x 10 12 cm -2 and Ni = 5 x 10 11 cm~ 2 (intermediate concentration regime), 
with a well defined hole impurity band and the valence band with its very small tail, not visible at this concentration. 
At right, (b) shows the SDF obtained for the two values of energy indicated on (a). The energy E\ lies in the middle 
of the impurity band, while E2 lies deep inside the conduction band. This is demonstrative of the fact that knowing 
the SDF at the Fermi level we can infer the character of our sample model as metallic or non-metallic according to 
the shape of the SDF. 

For each spin polarization the mobility edge is identifyed to the energy in which the SDF changes the shape from 
a sharp Lorentzian to a bell-shaped curve at k = 0. We have obtained this transition to be quite evident, the change 
occurring inside a narrow range of energy, as compared to the value of the Fermi level. Once the spin-dependent 
mobility edge is found for a given < M >, we obtain the density of extended carriers for each metallic sample: 



E 



E F 

dEM a {E). (14) 



The difference n s — n* represents the density of localized holes in the sample. We have, then, two conduction channels 
established: one, a "metallic" channel corresponding to the transport of the extended holes with energy lying between 
the mobility edge and the Fermi level, with density n*; another, that obtained by thermal activation of the localized 
states with activation energy A = \Ep — The total resistance is given by 

111 

7? — = r + r ' (15) 

-'•'total -^hopping -^metallic 

The first term on the right hand represents the resistance due to a hopping process, modeled by 

Chopping = -R00 exp(A/T). (16) 
The method we used to calculate the "metallic" resistivity is described in the next Section. 

III. RESISTIVITY OF THE "METALLIC" CHANNEL 

We assume quasi-two-dimensional carriers of density n* (arial density) inside a layered semiconductor structure 
were impurities (scatterer centers) are located inside specific layers with a density rii. The calculation is performed 
according to the Lei and Ting's 2 ^ force and energy balance equation formalism which has been used successfully in 
several heterostructures of different dimensionality^— both in the linear and in the non-linear regime. An advantage 
of this formalism is to provide a direct calculation of the self-averaged resistivity, instead of the conductivity. 

In the present case, a confining potential U(z) due to the structure interfaces breaks the symmetry in the z- 
direction, and in the absence of impurities and of an electric field, the single-particle states in the envelope function 
approximation are solutions of the Schrodinger equation: 

#o*™k(r, z) = £ nk #„k(r, z), (17) 

resulting into a free motion in the {x,y) plane and the envelope function </>„(z) in the z-direction: 

* nk (r,z) = -L e *"-0„(z), (18) 
VA 



with sub-band energies: 



h 2 k 2 

£nk = £n + i r ^- (19) 
2m* 



Equation p7|) is solved in the self-consistent Hartree-Fock approximation. 

Including scattering by impurities and the external electric field, the total Hamiltonian in the (x, y) plane can be 
separated into center of mass (CM) and relative coordinates of the N carriers, 

H = H C ra + H e + H e i, (20) 
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where 



H cm = K ^--NeK.E (21) 



He — ^nkC mk(T C m k(T (22) 

m,kcr 

H ei = U mn (q,z a )e-^< R+ ^cl k+(la c n]la (23) 



Vectors P and R. represent the momentum and in-plane position of the CM. We have used the creation and destruction 
operators c|„ kcr (c m k(T), and t/ m „(q, z a ) represents the potential matrix element of the interaction of an ion at position 
(r a ,z a ). It is given by: 

7p 2 f 

U mn (<l,z a ) = 2^ J dz<f>* ma (z)cj> n „(z)e-^-^, (24) 

where k is the GaAs dielectric constant. 

The configurational averaged drag force due to impurities depends on the polarization function: 

n mnCT;roW (k,k',q;t- t') = T([4 k ,_ qCT (t)c„ k v(i).4'k +qCT '(^)cn'k CT '(i')])o. (25) 
where T is the time ordering operator. Defining io d = q • v^, the averaged drag force is written as: 

fi(v d ) = X! q / dz a U mn (q,z a )U m > n >(-c[,z a )h(z a )^^n mna , m , n , a ,(k 7 k' ,q;w d ), 

mnt^'m'n' k<x k'u' 

with f(z a ) equals one if z a lies inside the impurity layer, zero otherwise, and 

/+oo 
dt e d ^ •'n nm(T!m ' n ' ( j' (k, k , q; t t) 
-oo 

We have neglected hole-hole correlation, therefore the polarization function is diagonal. In this case, 

ii(vd) = j- ^2\U mn (q)\ 2 qn^\m,n,(i,uj d ), (26) 

mnq 

with the effective potential due to impurities 

\Umn{q)\ 2 = J dz a \U mn {(\,z a )\ 2 h(z a ) (27) 



and 



(m, n, q, w d ) = 3 Y /o(^mk + q)-/o(U) (28) 

^ (mk+q — £nk — LOd — iO 

is the imaginary part of the non-interacting polarization function. In the stationary regime the balance of force on 
the system gives 

n* s eE + fi(v d )=0. (29) 

In that case, we have for the resistivity 



Assuming = VdX and isotropy, the resistivity becomes, using cod = VdQ x 

2 

y — \U mn {q)\ 2 n ( °\m,n,q,uj d ). (31) 



2hnfe 2 ^ uj d 1 

s mnq 
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In the Ohmic regime we take ujd — > 0, and 



Pi 



2hnfe 2 ^ ui d 

a mnq 



2 

J2 — I U ™ (l) 1 2 K n, q, 0) + Ilf } (m, n, q, 0)co d ] , 



(32) 



where 



Mo), n \ r ^n^ } (to, n, q, w d ) . 
(m,n, q,0) = — - — 



(33) 



IV. IMPURITY BAND AND TRANSPORT 



In what follows we present the calculation for the resistance of a double layer system, composed by two GaMnAs 
layers of width 20A immersed in a wide GaAs host. These two layers are separated by a GaAs spacer of width 
20A. First we assume that we have only the "metallic" channel, and we fix the carrier density while changing the 
magnetization. In that case, for a given carrier concentration n s , we observe a redistribution of the charge density 
in the z-direction as the average magnetization < M > decreases from < M >= 1 to < M >= 0, shown in the two 
extreme cases by Fig. |3l where m = 1. x 10 20 cm -3 and n s = 2x 10 13 cm~ 2 . The resistivity is expected to decrease 
as the average magnetization decreases, a consequence of the carriers occupying a wider region than that occupied by 
the impurities. In other words, we should expect the scattering cross-section to decrease towards the paramagnetic 
phase. 

At the same time as the sub-band occupation changes with (M), the Fermi level moves, as shown in the inset of Fig. 
Ujumping suddenly from a weakly occupied sub-band to a strongly occupied one. A similar behavior is observed as an 
effect of changing an applied magnetic field, so we should obtain the equivalent to the Shubnikov-de Haas oscillations. 
This is exactly what is observed in Fig. [4] showing the resistivity (notice the linear scale) as a function of the average 
magnetization. However, as (M) — ¥ several sub-bands become occupied, the spin-polarization disappears, several 
sub-bands of both spin orientation tends to contribute to the polarization function in equal foot. This provides an 
increase of the resistivity near (M) — 0. In conclusion, if we maintain the carriers density constant as we decrease 
the average magnetization, we obtain in our calculation oscillations of the resistivity in the range of strong magnetic 
order, followed by a minimum of the resistivity, and a hump as we approach (M) = 0. 

Next we consider the existence of an impurity band, appearing as a wide region where a non-zero density of states 
occurs below the edge of the valence band, superimposed to a tail of the valence band density of states. The signature 
of this impurity band shall be expressed by the existence of two transport channels, one due to eventually localized 
states, other to eventually extended states. The total density of particles, in our case holes, is assumed fixed. So, 
as the magnetization changes we have n*({M)) for the arial density of the extended states, and n s — n* for the 
localized states. The ratio n*/n s is obtained by the following ansatz. For a given pair of parameters (rii, n s ) labeling 
the sample, the spectral density function is used to determine the effective mobility edge according to the criteria 
described in Section [XT] Knowing the mobility edge we calculate directly the n*/n s ratio for a given m, as a function 
of the magnetization, as shown in Fig. [5] for three different samples. In sample #1 we have n s = 3. x 10 19 cm -3 and 
rii — 3.75 x 10 19 cm~ 3 , i.e, n s and n, are chosen very close to each other. Sample #2 has n s = 3.25 x 10 19 cm~ 3 and 
rii = 4.5 x 10 19 cm~ 3 . Finally sample #3 has n s = 3.25 x 10 19 cm~ 3 and rii = 7.5 x I0 19 cm~ 3 , n s as a small fraction 

Of Hi. 

Notice that in sample ^3 the ratio goes to approximately zero when the magnetization goes to zero, describing a 
sample with a very high resistance in the "metallic" channel (almost as an insulator) near Tc- The other two samples 
also show a decrease of the ratio as < M >—¥ 0, but remaining "metallic" near Tc. 

Now, given a sample n, we assume the ratio n*/n s to follow the same dependence as obtained in the 2D limit, and 
we take also the same activation energy, A = \Ep — fj,\, which is a reasonable approximation if we are not looking for 
quantitative comparison. On the other hand, in Eq. (|16[) is calculated by fitting the resistivity of samples #1, 
#2 and #5 in Ref. d to the function: 

R °° = ^78' ( 34 ) 

for a given Mn concentration x. We expected the dependence to be nearly proportional to the inverse of the average 
distance between Mn atoms. In fact we obtained (3 — 1.4. Then we interpolated the hopping resistance for the 
necessary Mn concentration. 

The calculation of the total resistivity for the same three samples is shown in Fig. [51 in logarithmic scale. Notice 
that curves for samples #1 and #2 are similar, showing a maximum of the resistivity at Tc, and a lower value of the 
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resistivity across the whole temperature range, as compared to sample #3, the one which behaves like an insulator 
near Tq. This latter shows a pronounced hump just before Tq and a flat region at low temperatures. This sample 
contains the clear signature of the impurity band. In fact, at low temperature, on the flat region, it shows a typical 
metallic behavior. Notice that the dependence of n* on < M >, together with the effect of finite temperature, fades 
the Shubnikov-de Hass oscillations. As we approach Tc the chemical potential gets close to the mobility edge, raising 
in importance the transport via thermal excitation of the localized states. 

To end our discussion about the character of metallic or non-metallic Gai-^Mn^As samples , Figs. OH] show 
schematically the relative positions of the impurity and valence band for three values of magnetization and three 
different situations as regard the classification of samples as metallic or nonmetallic. In this figure, the solid line 
curves represent, respectively, the DOS of anti-aligned and aligned spins with the magnetization. The vertical lines 
represent their mobility edges and the Fermi level, which moves to the left as the average magnetization goes to zero, 
once we fixed the left edge of the anti-aligned DOS. 

In Fig. [7] we take a sample in which the Fermi level lies to the left of the anti-aligned mobility edge. This samples 
is an insulator even at < M >= 1. As < M > decreases the splitting decreases and the Fermi level goes to the left, 
keeping all the states localized. In Fig. [8] the carrier density is high enough so that even with the Fermi level going 
to the left as the magnetization decreases, it is maintained to the right of the anti-aligned mobility edge. This sample 
behaves as a metal even at < M >= 0, like samples # I and # 2. Fig. [9] represents a sample which is "metallic" 
at < M >= 1, but as the magnetization decreases the Fermi level moves to the left of the spin anti-aligned mobility 
edge, and the sample behaves as non-metallic below a certain < M > . 

It is remarkable that the schematic case of Fig. ® is typically the situation observed by the calculation of the 
resistivity of sample # 3 in Fig. [51 containing the signature of the impurity band. Notice that, since we are treating a 
quasi-two-dimensional sample, the density of states just after the tail of the valence band has a very sharp derivative, 
as can be seen in the real calculation of Fig. [5J This means that the Fermi level is expected to be located just after 
this tail, if the sample has a metallic behavior. In that case, depending on the carrier concentration, when the Fermi 
level moves to the left as the magnetization decreases, it may occupy the region below the valence band tail, in other 
words, a position already inside the impurity band. In that case, a change of sign on the derivative of the DOS should 
be observed. It is important to mention that this is in line with the explanation for the observation of anomalous 
Hall effects performed in GaMnSb in Ref. |33|. 



V. CONCLUSION 



These calculations present theoretical evidence of the existence of an impurity band contribution to the transport 
when the curve of the resistivity versus temperature has a typical metallic behavior at low temperature and a pronoun- 
ced hump with a maximum just before the transition temperature. Although the calculations have been performed in 
a double layer structure, where the effects of spin-polarization are enhanced, and, in consequence, the occurrence of 
metallic to a non-metallic behavior as the critical temperature is approached, we believe that our conclusion can be 
extended from the digital up to the epilayer case, and is an additional evidence, together with optical measurements, 
of the importance of a hole impurity band in GaMnAs. It is worthwhile to add that the existence of an impurity 
band creates not only two transport channel, but two interaction mechanism between localized magnetic moments 
as well: one via localized states (possibly a magneto-polaron mechanism), another via extended states (RKKY, for 
instance). Samples like the one represented schematically by Fig. [5] would have both mechanism acting up to a certain 
value of the temperature, above which the low average magnetization makes the mobility edge to lie above the Fermi 
level, and the "extended" states channel ceases acting. From that point on, the magnetization decreases faster with 
temperature. We believe this is the reason for which it is hard to approximate the curve magnetization as a function 
of temperature in layered GaMnAs structures by a simple Brillouin function^ 
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Figura 1: Self-energy diagrams, within the multiple scattering approximation, and the Dyson equation for the averaged GF. 
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Figura 2: Typical spectral density function for localized states (like that with energy E\) and extended states (like that with 
energy E2). 
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Figura 3: (color online) Charge distribution across the z-direction for (M) = and (M) = 1. 
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Figura 4: (color online) Resistivity as a function of the average magnetization for fixed carrier density. The scale is chosen in 
such a way that the temperature corresponding to a given magnetization grows to the right. The inset shows the index of the 
sub-band containing the Fermi level. 
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Figura 5: (color online) Dependence of the ratio n* a /n a with < M > for three different samples; sample #1: n s = 3. x 10 19 cm 3 
and n x = 3.75 x 10 19 cm _3 ; sample #2: ?i s = 3.25 x 10 19 cm -3 and m = 4.5 x 10 19 cm _3 ; sample #3: n a = 3.25 x 10 19 cm -3 and 
m = 7.5 x 10 19 cm -3 . 
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Figura 7: (color online) Schematic classification of samples as metallic or nonmetallic. The downward (upward) arrow indicates 
the anti-aligned (aligned) spin impurity band and its corresponding mobility edge. The Fermi levels lies, in this case, to the 
left of the anti-aligned mobility edge and the sample is non-metallic in the whole range of average magnetization. 
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Figura 8: (color online) As before, but the Fermi level is kept to the right of the anti-aligned mobility edge, even as < M > 
goes to zero. This is the case of a sample which is metallic for the whole range of the magnetization. 




Figura 9: (color online) As before, but in this case the mobility edge lies to the right of the anti-aligned mobility edge 
< M >= 1, but crosses to its left as M goes to zero. 



